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ABSTRACT 

We present the first simultaneous, multi-wavelength observations of an L dwarf, the L3.5 candidate brown 
dwarf 2MASS J00361617+1821 104, conducted with the Very Large Array, the Chandra X-ray Observatory, and 
the Kitt Peak 4-m telescope. We detect strongly variable and periodic radio emission (P = 3 hr) with a fraction 
of about 60% circular polarization. No X-ray emission is detected to a limit of Lx/^boi ^2x 10~ 5 , several 
hundred times below the saturation level observed in early M dwarfs. Similarly, we do not detect Ha emission 
to a limit of L// a /Lb i ^2x 10~ 7 , the deepest for any L dwarf observed to date. The ratio of radio to X-ray 
luminosity is at least four orders of magnitude in excess of that observed in a wide range of active stars (including 
M dwarfs) providing the first direct confirmation that late-M and L dwarfs violate the radio/X-ray correlation. 
The radio emission is due to gyrosynchrotron radiation in a large-scale magnetic field of about 175 G, which is 
maintained on timescales longer than three years. The detected 3-hour period may be due to (i) the orbital motion 
of a companion at a separation of about five stellar radii, similar to the configuration of RS CVn systems, (ii) 
an equatorial rotation velocity of about 37 km s" 1 and an anchored, long-lived magnetic field, or (iii) periodic 
release of magnetic stresses in the form of weak flares. In the case of orbital motion, the magnetic activity may 
be induced by the companion, possibly explaining the unusual pattern of activity and the long-lived signal. We 
conclude that fully convective stars can maintain a large-scale and stable magnetic field, but the lack of X-ray and 
Ha emission indicates that the atmospheric conditions are markedly different than in early-type stars and even M 
dwarfs. Similar observations are therefore invaluable for probing both the internal and external structure of low 
mass stars and sub-stellar objects, and for providing constraints on dynamo models. 

Subject headings: stars: low mass,brown dwarf-stars: activity-stars: radio emission-stars: magnetic 
fields-radiation mechanisms: nonthermal 
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spectral type M8 (e.g., Gizis et al. 2000; West et al. 2004). 
Moreover, in both bands the level of activity increases with 
rotation velocity, up to a saturation level, at v ~ 5 km s" 1 , of 
Lx/L ho \ « 10" 3 (Rosner, Golub & Vaiana 1985; Fleming et al. 
1993; Pizzolato et al. 2003) and L Ha /L hoi « 10" 3 5 (Mohanty 
et al. 2002). These observations have been interpreted in the 
context of coronal and chromospheric heating by dissipation 
of magnetic fields, which in turn are powered by an internal dy- 
namo. In stars earlier than M3 the dynamo is thought to be pow- 
ered by shearing motions at the radiative-convective transition 
zone (the so-called afl dynamo). The increase in activity with 
decreased mass and higher rotation velocity is understood as a 
dependence of the dynamo on the Rossby number, Ro = P/t c , 
where P is the rotation period and t c is the convective turnover 
time (e.g., (Soderblom et al. 1993)). The saturation effect is 
still not fully understood. 

Surprisingly, the onset of full convection and the expected 
breakdown of the aQ dynamo at spectral type M3 is not ac- 
companied by obvious changes in the level of Ha and X-ray 
emission or the rotation-activity relation. This has been at- 
tributed to a growing contribution from an a 2 (Raedler et al. 
1990) or a turbulent dynamo (Durney, De Young & Roxburgh 
1993) in which the shearing needed to generate and amplify the 
fields comes from turbulent motion associated with the inter- 
nal convection. However, beyond spectral type M7 there is a 
precipitous drop in Ha and X-ray persistent activity, and only 
a handful of objects exhibit flares, <~ 7% in late-M dwarfs and 
~ 1% in L dwarfs (e.g., Reid et al. 1999; Gizis et al. 2000; 
Rutledge et al. 2000; Liebert et al. 2003; West et al. 2004 
and see Figure 14). Furthermore, the activity-rotation relation 
no longer holds in late-M and L dwarfs with many rapid rota- 
tors (v J> 20 km s" 1 ) exhibiting no discernible activity (Basri & 
Marcy 1995; Mohanty & Basri 2003). The decrease in activ- 
ity may be due to a quenching of the turbulent dynamo and a 
transition to small-scale and short-lived fields 17 . Moreover, the 
increasingly neutral atmospheres of late-M and L dwarfs may 
hamper the dissipation of the fields (Mohanty et al. 2002). In 
effect, the increased rate of collisions between charged and neu- 
tral particles may decouple the magnetic field from the atmo- 
sphere, thereby suppressing the build-up of magnetic stresses 
and the heating of the chromosphere and corona. 

The above discussion suggests that radio emission, produced 
by interaction of relativistic electrons with the magnetic field, 
should also be suppressed in late-M and L dwarfs. Observa- 
tionally, the empirical relation between the radio and X-ray lu- 
minosities of a wide range of coronally-active stars, including 
many M dwarfs (Guedel & Benz 1993; Benz & Guedel 1994), 
along with the few X-ray detections and upper limits obtained 
to date, suggest an expected flux level below 0.1 /xJy (i.e., sev- 
eral hundred times lower than the detection threshold of the 
Very Large Array). 

However, radio observations of late-M and L dwarfs present 
a decidedly different picture. Berger et al. (2001) (hereafter, 
B01) detected radio flares and persistent emission from the 
M9.5 brown dwarf LP944-20 four orders of magnitude brighter 
than expected based on X-ray observations (Rutledge et al. 
2000; Martin & Bouy 2002) and the radio/X-ray correlation. A 
subsequent survey of twelve additional objects (Berger 2002; 



hereafter, B02) revealed flares and persistent emission from 
three objects ranging from M8.5 to L3.5. In all cases the in- 
ferred magnetic fields strengths were of the order of 100 G. 
In addition, of the three detected objects, the two which have 
been observed in X-rays again violated the radio/X-ray corre- 
lation by orders of magnitude. Recently, Putman & Burgasser 
(2003) detected radio emission from four out of seven M7 to 
L4.5 dwarfs they observed with Australia Telescope Compact 
Array, of which one exhibited a strong flare. Thus, the inci- 
dence rate of radio emission in late-M and L dwarfs may be as 
high as 40%, much larger than in Ha or X-rays. 

Beyond the individual properties inferred from the unex- 
pected radio emission, two puzzling trends have emerged. First, 
the fraction of persistent radio emission compared to the bolo- 
metric luminosity appears to increase with later spectral type, 
the exact opposite of the trend observed in Ha and X-ray activ- 
ity. Second, the radio bright objects are those with high rotation 
velocities (vsim ;> 20 km s" 1 ), while significant non-detections 
are associated with slow rotators (vsin; <J 10 km s -1 ). This im- 
plies that the rotation-activity relation holds in the radio, despite 
being broken in Ha and X-rays. 

The glaring disparity between the trends observed in the var- 
ious activity indicators makes it clear that a complete picture 
of magnetic activity in late-M and L dwarfs requires simul- 
taneous, multi-wavelength observations. These observations 
can furthermore constrain turbulent dynamo models, which de- 
spite significant progress are still under-developed and under- 
constrained. Here we present the first such observations for any 
L dwarf, the L3.5 object 2MASS J00361617+1821 104 (here- 
after, 2M0036+18). In §2 we provide the details of our radio, 
X-ray, and Ha observations. We summarize the results in §3 
and address the radio emission mechanism and derive the mag- 
netic field properties in §4. The violation of the radio/X-ray 
correlation and its implications are discussed in §5. A possible 
periodicity observed in the radio emission is assessed in §6 and 
implications for the magnetic field generation mechanism are 
drawn in §7. Finally, we summarize the current trends in radio, 
Ha and X-ray emission from M and L dwarfs in §8. 

2. OBSERVATIONS 

We targeted the L3.5 dwarf 2M0036+18 (Reid et al. 1999) 
due to its relative vicinity and previous observations of radio 
activity. The object is located at a distance of 8.8 pc (Dahn 
et al. 2002), has a bolometric luminosity, Lboi ~ 10~ 3,97 L Q 
(Leggett et al. 2002; Vrba et al. 2004), a rotation velocity, 
vsin/ w 15 ±5 km s -1 , measured from high-resolution optical 
spectra (Schweitzer et al. 2001), an inferred temperature of 
about 1900 K (Vrba et al. 2004), an inferred radius of about 
0.09 ± 0.01 R (Dahn et al. 2002), a mass of about 0.076 M Q 
inferred from the surface gravity, logg « 5.4 (Schweitzer et al. 
2001), and an inferred age of at least 1 Gyr (e.g., Burrows et al. 
2001). Thus, 2M0036+18 is an object located near or below 
the hydrogen burning limit. 2M0036+18 has been detected in 
previous radio observations carried out in September and Oc- 
tober of 2001, with a strong ~ 20-min flare, variable persistent 
emission, and a high fraction of circular polarization (B02). 
The observations presented here were carried out simultane- 

17 Recent work by Dobler, Stix & Brandenburg (2004) on magnetic field generation in convective rotating stars suggests that large-scale fields can be generated and 
maintained on relatively long timescales. However, these authors use physical conditions relevant for an M5 dwarf (and with an unusually high surface temperature), 
while substantial changes in Ha and X-ray activity are manifested only for objects later than M7, and we are interested in mid-L dwarfs here. 
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ously in the radio, X-ray, and optical on September 28, 2002. 
Radio observations commenced at 01 :49 UT and ended at 09:39 
UT (28.2 ks), X-ray observations covered the range 01:10 to 
07:17 UT (22 ks), and optical spectroscopic observations cov- 
ered the range 02:47 to 07:31 UT (17 ks). 

2.1. Radio 

Very Large Array (VLA 18 ) observations were conducted si- 
multaneously at 4.86 and 8.46 GHz in the standard continuum 
mode with 2 x 50 MHz contiguous bands at each frequency. 
We used fourteen antennas at 4.86 GHz and thirteen antennas at 
8.46 GHz, with a staggered configuration in each of the three ar- 
ray arms. Eight-minute scans on 2M0036+18 were interleaved 
with 50 s scans on the phase calibrator J0042+233. The flux 
density scale was determined using the extragalactic source 3C 
48 (J0137+331). 

The data were reduced and analyzed using the Astronomical 
Image Processing System (AIPS; Fomalont 1981). The visi- 
bility data were inspected for quality, and noisy points were 
removed. To search for flares, we constructed light curves 
using the following method. We removed all the bright field 
sources using the AIPS/IMAGR routine to CLEAN the region 
around each source (with the exception of 2M0036+18), and 
the AIPS/UVSUB routine to subtract the resulting source mod- 
els from the visibility data. We then plotted the real part of the 
complex visibilities at the position of 2M0036+18 as a func- 
tion of time using the AIPS/UVPLT routine. The subtraction of 
field sources is necessary since their sidelobes and the change 
in the shape of the synthesized beam during the observation re- 
sult in flux variations over the map, which may contaminate 
any real variability or generate false variability. We repeated 
this procedure for the background sources as a check on the in- 
trinsic noise properties of the maps. The resulting light curves 
are shown in Figures 1 and 2. 

We conducted follow-up observations with the VLA on 
2005, Jan. 10 and 11 UT at 4.86 GHz. A total of ten hours 
were obtained using the same phase and flux calibrators as in 
the Sep. 2002 observation. We followed the same procedure 
outlined above to to produce light curves (Figures 8 and 9) 

2.2. X-rays 

The observations were made with the Chandra/ ACIS-S3 
(backside-illuminated chip). 2M0036+18 was offset from the 
on-axis focal point by 4' to mitigate pileup in the event of a 
flare. Data were analyzed using CIAO v3.1. We extracted data 
within a 4" circle centered on the source position, selecting 
counts in a low energy band (0.2-3 keV) and a total energy band 
(0.2-10 keV). Sources were found using celldetect. To es- 
timate the background we used annuli centered at the source 
position, excluding other point sources detected in the obser- 
vation. We adopt an energy conversion factor of 1 count = 
4.5 x 10" 12 erg cm -2 s" 1 (0.2-8 keV), consistent with a kT=l 
keV Raymond-Smith plasma. 

No X-ray source is detected at the optical position of 
2M0036+18. The number of counts in the low and total en- 
ergy bands were 3 and 5, respectively, while we expected 2 and 
4 counts from the background. There is no evidence for flaring 
activity, or deviation from the Poisson count rate distribution 

18 The VLA is operated by the National Radio Astronomy Observatory, a fi 
Associated Universities, Inc. 



expected from a constant background. For example, the short- 
est time-span which includes 3 counts is 6.4 ks. For a detection 
of 5 counts, the 90% confidence limit (found from an average 
count rate of <9 counts per 21 .5 ks) on the time-averaged flux is 
< 1 x 10~ 15 erg cm" 2 s" 1 . We observed no more than 1 count per 
1 hr time period, corresponding to a 90% confidence limit on a 
1 hr average peak flux (<3 counts/hr) of < 4 x 10~ 15 erg cm -2 
s" 1 . At the distance of 2M0036+18 the upper limit on persistent 
emission is L x <; 9.3 x 10 24 erg s" 1 , or log(Lx/£-boi) <, -4.65. 

2.3. Optical 

2M0036+18 was monitored spectroscopically using the 
Multi-Aperture Red Spectrometer (MARS) mounted on the 
Mayall 4-meter telescope at Kitt Peak National Observatory. A 
series of forty-eight 300-s exposures covering 6400- 10, 100A 
with a resolution of 3A pix" 1 were obtained in clear and photo- 
metric conditions with reasonable seeing (<~1.2", although the 
initial observations were made at airmass exceeding 2). All 
of the spectra were flat-fielded, extracted and wavelength cal- 
ibrated using standard IRAF routines. Observations of Feige 
110 were used to set the flux scale. 

None of the spectra show evidence for Ha emission. The 
individual 300-s spectra have signal-to-noise of about 10 to 
15 at these wavelengths, leading to an upper limit on the Ha 
equivalent width of w 1.5A (or F Ha <; 5 x 10" 17 erg cm" 2 s" 1 ). 
Combining all the individual spectra (Figure 3), we find no ev- 
idence for emission exceeding an equivalent width of 0.4A, or 
F Ha & 1 x 10" 17 erg cm" 2 s _1 . At the distance of 2M0036+18 
this translates to L Ha & 9.3 x 10 22 erg s"\ or log(L Ha /L bol ) £ 
-6.65. 

Separately, on 2004 December 5 we obtained two 10-min 
high-resolution spectra with the Magellan Inamori Kyocera 
Echelle (MIKE) Spectrograph, separated by 103 min. The pur- 
pose of this observation was to measure a possible radial veloc- 
ity induced by a close-in companion (see §6). A 0.35" slit was 
used with 2x2 binning, providing a spectral resolution of about 
27,000 in the red or about 0.11 A per pixel. The spectra were 
flat-fielded, rectified, wavelength calibrated and extracted using 
a MIKE data reduction pipeline. The rms wavelength residuals 
were of the order of 0.005A and the signal-to-noise in the rel- 
evant echelle orders was approximately ten and six in the first 
and second spectrum, respectively. 

3. BASIC PROPERTIES OF THE RADIO EMISSION 

While neither X-ray nor Ha emission were detected from 
2M0036+18, radio emission was detected with an average flux 
of 259 ± 19 fiJy and 134 ± 16 ^Jy at 4.9 and 8.5 GHz, re- 
spectively (Figure 1). The resulting spectral index (F v oc v a ) 
is a = -1.2 ±0.3, typical of radio emission observed from M 
dwarfs (Gudel et al. 1993). The emission is strongly circu- 
larly polarized with an average of f c w -73 ± 8% at 4.9 GHz 
and f c « -60 ± 15% at 8.5 GHz. The negative sign indicates 
that the polarization direction is left-handed. We place a limit 
of <J 20% (3cr) on the fraction of linear polarization at 4.9 GHz 
and <J 35% at 8.5 GHz. In the follow-up observations we de- 
tected emission from 2M0036+18 at a level of 152 ±9 /uJy, with 
a fraction of circular polarization, f c w -46 ± 7%. 

For comparison, two past observations of 2M0036+18 at 8.5 

of the National Science Foundation operated under cooperative agreement by 
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GHz (B02) revealed a similar flux level, 135 and 330 /iJy, but 
a lower fraction of circular polarization, -35% and -13%. One 
of the two observations also uncovered a 20-min flare with a 
peak flux of about 720 yuJy and a fraction of circular polariza- 
tion, f c sa -62%. No such flares are observed here suggesting 
that their rate of occurrence, taking into account all three obser- 
vations, is < 0.04 hr" 1 . The flare duty cycle, defined as the flare 
duration relative to the total observing time, is <J 1.5%. 

The brightness temperature of the radio emission, an impor- 
tant discriminant of the emission mechanism, is given by 

T b « 2 x 10 9 F V:m]y ^ Hz d 2 pc (R/R,r 2 K, (1) 

where R s /Rj w 1 is the radius in units of Jupiter radii. We find 
that Ti, ~ 2 x 10 8 K if the emission is produced in the corona 
at a radius of R ~ 3R S ~ 3Rj (Linsky & Gary 1983; Burrows, 
Hubbard & Lunine 1989; Leto et al. 2000) and with a covering 
fraction of 100%. On the other hand, if the emission is con- 
fined to a significantly smaller region (e.g., a coronal loop with 
R <J Q.IRs) then 7], may exceed 10 11 K. In the case of gyro- 
magnetic radiation (see §4.1), the inverse Compton catastrophe 
limit, Tb <J 10 12 K, defines a minimum size for the emission 
region of about 0.08/?, w 5.0 x 10 8 cm. 

4. EMISSION MECHANISM 

Several radiation mechanisms can give rise to the observed 
radio emission (e.g., Giidel 2002): Bremsstrahlung radia- 
tion, gyromagnetic radiation, and coherent processes (plasma 
radiation or electron cyclotron maser). We can rule out 
bremsstrahlung radiation in this case since the expected spec- 
trum is flat in the optically thin regime, whereas here a w —1.2. 

Coherent emission processes, which have been invoked in 
cases of high brightness temperature (7], > 10 12 K) solar and 
stellar flares with strong circular polarization, are also unlikely 
in this case. This is primarily because the emission is expected 
to have a narrow bandwidth around the fundamental or second 
harmonic of the plasma frequency, v p « 9000«y 2 Hz, or the cy- 
clotron frequency, v c k> 2.8 x 10 6 Z? Hz, whereas here emission 
is detected at both 4.9 and 8.5 GHz. In addition, the typical 
emission timescale of coherent radiation flares is much shorter 
than the 7.8-hour duration of the radio emission detected in this 
case. 

4. 1 . Gyromagnetic Radiation 

Gyromagnetic radiation can arise from a thermal or non- 
thermal distribution of electrons, with the latter typically as- 
sumed to follow a power law, n(j) oc j~ p , above some cut- 
off, 7 m . The emission properties further depend on the har- 
monic number, s, or equivalently the typical energy of the elec- 
trons 7 = s 1 / 3 . Cyclotron emission is produced when s <J 10 
(non-relativistic electrons, typically thermal distribution), gy- 
rosynchrotron when s~ 10- 100 (mildly relativistic electrons), 
and synchrotron when 100 (relativistic electrons). The fre- 
quency at which the bulk of the radiation is emitted is given by 
v m « 2.8 x 10 6 sB Hz, where B is the magnetic field strength. 

Cyclotron radiation is unlikely in this case since the emis- 
sion is expected to be concentrated in narrow bands at the har- 
monic frequencies, as opposed to the detected signal at 4.9 and 
8.5 GHz. Similarly, gyrosynchrotron emission from a thermal 
plasma is unlikely since the expected optically thin spectral in- 
dex, a w -8, is significantly steeper than the observed value. 



Finally, synchrotron emission is unlikely because the expected 
degree of circular polarization is strongly suppressed compared 
to the value of about 70% measured here. Therefore, the most 
likely emission mechanism in the case of 2M0036+18 is gy- 
rosynchrotron radiation from a non-thermal distribution of elec- 
trons. 

The emission properties in the gyrosynchrotron case depend 
both on the fundamental source parameters - magnetic field B, 
electron density n e , and line-of-sight source scale L - and the 
angle between the line of sight and the magnetic field, (Dulk 
& Marsh 1982). Since the latter cannot be inferred from the 
present data we take 9 = n/3 as a representative value. The 
spectral index in the optically thin case, a w 1.2-0.9/?, indi- 
cates that p « 2.7. This value is similar to those found in sev- 
eral early M dwarfs (Gudel et al. 1993). Given the value of 
p, the fraction of circular polarization at 4.9 GHz is given by 
f c « 0.025B 53 « 0.6, from which we infer B « 400 G. For 
the full range of possible values the magnetic field strength 
ranges from about 80 to 1 .9 x 10 3 G. These values indicate that 
the harmonic number is s < 22, with the gyrosynchrotron limit, 
s > 10, corresponding to B < 175 G, or 6 < ir/4. Below we use 
the latter value for the magnetic field strength. 

The frequency at which the gyrosynchrotron spectrum peaks 
is given by v m w 6.1 x 10 5 (n e L) 24 < 4.86 x 10 9 Hz, while the 
flux in the optically thin regime is F v sa 1 . 1 x 10~ M n e LR 2 « 205 
fijy, using the average flux at 4.9 GHz. Solving for the un- 
known column density and radius of the emission region we 
find n e L < 1 .8 x 10 16 cm -2 and R « 1 x 10 10 cm, or about 1 AR S . 
The resulting brightness temperature is Tb w 7 x 10 8 , consis- 
tent with the interpretation of the radio emission as gyrosyn- 
chrotron radiation. Assuming that L = R, the electron density is 
n e « 1.6 x 10 6 cm" 3 . 

Several conclusions can be drawn from this discussion. First, 
the magnetic field strength is larger than the average solar field 
or the magnetic field of Jupiter, but is somewhat lower than the 
~ 1 kG fields inferred in some active early M dwarfs (Smith 
et al. 1975; Saar & Linsky 1985; Haisch, Strong & Rodono 
1991; Johns-Krull & Valenti 1996; Stepanov et al. 2001). We 
note however, that the surface magnetic field of 2M0036+18 
may be an order of magnitude larger than that inferred for the 
emission region if L w R « IAR S , placing it in the same range 
of early M dwarfs. Second, the radio emission, and in particular 
the high fraction of circular polarization, requires a large-scale, 
long-lived ordered magnetic field with a large covering frac- 
tion. This is again similar to the conditions inferred for early M 
dwarfs (Saar & Linsky 1985; Johns-Krull & Valenti 1996). Fi- 
nally, the density of emitting electrons is similar to that inferred 
in the emission zones (coronae) of early M dwarfs (Stepanov 
et al. 2001; Osten et al. 2004). Thus, the physical conditions 
in the L3.5 2M0036+ 18 resemble those found in active early M 
dwarf, but the Ha and X-ray emission are strongly suppressed. 

We finally note that the plasma-beta, = 2fioP/B 2 , of the 
corona, using n ~ 1 .6 x 10 6 cm" 3 , a typical temperature r~2x 
10 6 K (Giampapa et al. 1996), and B ~ 175 G, is ~ 3 x 10" 8 . 
The typical photospheric pressure of an L dwarf is P ~ 3 x 10 5 
dyne cm" 2 (Burrows et al. 2001), resulting in ~ 2. These con- 
ditions are again similar to those found in M dwarfs (Giampapa 
et al. 1996) and the Sun. Thus, even small perturbations in 
the coronal field of 2M0036+18 are capable of supporting the 
pressure gradients produced by the radio-emitting material. 
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5. VIOLATION OF THE RADIO/X-RAY CORRELATION 

The simultaneous radio and X-ray observations allow us for 
the first time to directly investigate the radio/X-ray correla- 
tion in an L dwarf. Similar observations of coronally active 
stars (including the Sun) reveal a tight correlation between 
L R and L x (Figure 4; Guedel & Benz 1993; Benz & Guedel 
1994). This correlation holds for both flaring and persistent 
emission, and for single stars and binary systems. For M 
dwarfs the correlation holds to spectral type M7 (Gudel et al. 
1993). The persistent emission follows L R w 3 x 10 _16 L X Hz -1 
(Guedel & Benz 1993), and extends over six orders of magni- 
tude in L R . The relation for radio and X-ray flares is non-linear, 
L* oc 1.3 x 10- 26 L° 73 Hz" 1 (Benz & Guedel 1994), and holds 
over eight orders of magnitude in L R . 

For 2M0036+18 we find L R w 2.5 x 10 13 erg s" 1 Hz" 1 and 
L x < 9.7 x 10 24 erg s _1 . Thus, L R /L X £ 2.6 x 1(T 12 Hz" 1 is at 
least 8600 times higher than predicted by the radio/X-ray corre- 
lation (Figure 4). A similar result, from non-simultaneous data, 
was obtained for the M9.5 brown dwarf LP944-20 (B01) with a 
ratio that exceeded the expected value by a factor of 10 4 (flares) 
and > 8500 (persistent emission). Similarly, for the M9.5 dwarf 
BRI 0021-0214 the radio luminosity exceeded the X-ray upper 
limit by at least a factor of 10 3 (B02). Thus, the radio and X-ray 
luminosities appear to de-correlate over a narrow range in spec- 
tral type (between about M7 and M9) and remain uncorrected 
at least to mid-L (Figure 4). Incidentally, this is the same range 
over which both the fraction of Ha and X-ray active sources 
and the strength of the activity drop precipitously. 

The violation of the radio/X-ray correlation is surprising 
given that the magnetic and emission region properties derived 
for 2M0036+18 are not dissimilar from those found on early 
M dwarfs (§4.1). This suggests that either the origin of the 
magnetic field is different than in M dwarfs or the atmospheric 
conditions are not conducive for the support of a chromosphere 
and/or a corona. 

The radio/X-ray correlation has been interpreted in the con- 
text of magnetic reconnection. In this scenario, the flares are 
produced when coronal magnetic loops reconnect, release en- 
ergy, and create a current sheet along which ambient electrons 
are accelerated. The accelerated electrons drive an outflow of 
hot plasma into the corona as they interact with and heat the 
underlying chromospheric material. The interaction of the out- 
flowing plasma with the electrons produces X-ray emission via 
the bremsstrahlung process (Neupert 1 968 ; Hawley et al. 1 995 ; 
Guedel et al. 1996). This so-called Neupert effect points to a 
causal connection between particle acceleration, which is the 
source of radio emission, and plasma heating, which results in 
X-ray emission. Thus, the X-ray thermal energy should sim- 
ply be related by a constant of proportionality to the integrated 
radio flux. 

Guedel & Benz (1993) discuss this correlation in terms of the 
relation L R /L X « lO~ 22 B 2 - 5 (a/b)T () (l + a) Hz" 1 , where a is the 
fractional efficiency of accelerating electrons, b is the fraction 
of coronal energy radiated in X-rays, tq is the average lifetime 
of the low-energy electrons, and a is the power law relating the 
electron lifetime to its energy. The observed tight correlation 
suggests that the combination B 2 5 (a/b)TQ is nearly unchanged 
for a wide range of stars. Since in the case of 2M0036+18 and 
the other radio active late-M and L dwarfs L R /L X is higher by 
several orders of magnitude while the inferred value of B is 



not unusual compared to M dwarfs, this indicates that either 
the fraction of energy emitted in X-rays, b, is suppressed by a 
factor of ~ 10 4 , or the typical lifetime of the electrons, to, is 
longer by a similar factor. In particular, for a/b of order unity 
and 8 s 175 G, we find a mean electron lifetime t ~ 1 d. In this 
case, the X-ray emission is suppressed because the electrons do 
not lose their energy to coronal heating on a sufficiently rapid 
timescale. 

On the other hand, the violation of the radio/X-ray correla- 
tion may be due to an inefficient production of X-rays, i.e., a 
low value of b. In the case of L dwarfs this may be explained 
by the increasingly neutral atmospheres and the reduction in 
the density of ions available for X-ray emission (Mohanty et 
al. 2002). However, it is not clear if this can account for the 
late-M dwarfs LP944-20 and BRI 0021-0214, or for the sud- 
den transition between spectral types M7 and M9. More likely, 
the reduction in X-ray luminosity, or the over-production of ra- 
dio emission, may be related to the magnetic structure across 
the atmosphere, chromosphere and corona or efficient trapping 
of the electrons (i.e., long lifetime). We note that recent radio 
and X-ray observations of the dM4.5e star EV Lac also reveal 
a breakdown of the Neupert effect, which Osten et al. (2004) 
interpret in the same context discussed here, namely efficient 
trapping or a low value of b. 

6. VARIABILITY AND PERIODICITY 

To this point we have discussed the properties of the inte- 
grated radio emission during the entire 7.8 hour simultaneous 
observation and the follow-up observations. However, from 
Figures 1, 2, 8, and 9 it is clear that both the total intensity 
and circularly polarized emission at 4.9 GHz are strongly vari- 
able. For the total intensity light curve the reduced \ 2 assum- 
ing a constant flux is w 2.5-7.5 for a time resolution ranging 
from 5 min to 1 hour; for the circular polarization light curve 
X 2 w 1.8-2.8 over the same range of time resolutions. It is 
not clear if the variability of the weaker emission at 8.5 GHz is 
significant, with \ 2 r ~ 1 — 1.7 depending on the time resolution. 

We are confident that this variability is not an observational 
artifact for two reasons. First, the light curves of two field 
sources, constructed in the same manner as that of 2M0036+18 
(see §2.1), do not exhibit significant variability, with % 2 w 1.2 
and « 1.8 for a time resolution of 10 min. In fact, the re- 
duced \ 2 values for the field sources are over-estimated com- 
pared to those for 2M0036+18 since, unlike 2M0036+18, these 
sources are spatially extended and positioned away from the 
phase center. As a result, changes in the synthesized beam dur- 
ing the observation affect their light curves more significantly 
than that of 2M0036+ 1 8 . Second, the variability of 2M0036+ 1 8 
is significant in circular polarization where no field sources, 
which can introduce false variability, produce detectable emis- 
sion. It is important to note that in the two past observations of 
2M0036+18, lasting 150 and 180 minutes, the persistent emis- 
sion was also variable (B02). 

In addition to the variability, it appears from Figures 2, 8, 
and 9 that the radio emission at 4.9 GHz is periodic. The auto- 
correlation function (ACF) of both the total intensity and circu- 
lar polarization at a wide range of time resolutions shows a clear 
peak at 3 hours (Figures 5 and 10). The ACFs of the two field 
sources described above are flat, as expected for non-periodic 
sources. 

We assess the statistical significance of the 3-hour periodic- 
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ity using two other methods which provide a measure of the 
power spectrum and are therefore more easily calibrated: The 
Lomb-Scargle (LS) periodogram (e.g., Press et al. 1992), and a 
one-dimensional CLEAN algorithm (Roberts, Lehar & Dreher 
1987). The latter method has the added advantage that artifacts 
and spectral leakage which arise from the sampling function 
and the finite length of the observation are removed. 

We assess the significance of the results using the Monte 
Carlo method. Specifically, we construct both randomized ver- 
sions of the 4.9 GHz light curves, and simulated light curves 
with a (x 2 ) variability similar to that of 2M0036+18. Depend- 
ing on the time resolution we use 50,000 to 100,000 light 
curves in each set and perform the same procedures as on the 
real data. The significance of peaks in the LS and CLEAN 
power spectra of the real light curve are determined as the frac- 
tion of simulated light curves which produce stronger peaks. 
We find that both the randomized and simulated sets provide a 
similar measure of the significance, as do the LS and CLEAN 
procedures. A representative LS periodogram is shown in Fig- 
ure 6. 

The significance of peaks in the LS and CLEAN power spec- 
tra does depend on the choice of time resolution. This is ex- 
pected since as the period is sampled more coarsely the strength 
of the signal is expected to decrease. We use time resolu- 
tions ranging 19 from 2 to 30 min (Figure 7), and find that the 
strongest peak always occurs at a frequency of 0.33 hr" 1 , or 
a period of 3 hours. For the simulated light curves the peaks 
that exceed the value for 2M0036+18 occur over the entire fre- 
quency range. We also find that the significance of the peri- 
odicity in the circular polarization light curve is consistently 
higher than for the total intensity light curve. This is expected 
because the contamination from field sources is non-existent in 
circular polarization. For the range of time resolutions investi- 
gated we find that the significance of the 3-hour period peaks 
at w 99.999% for St = 2 min. This result is confirmed by the 
follow-up observations, and a phased light curve of the data 
from 2005, Jan. 10 and 1 1 UT, folded with a period of 184 min 
is shown in Figure 1 1 . 

Thus, the radio emission from 2M0036+18 is periodic with 
P = 3 hours. This period is maintained on timescales of at least 
two years. If related to the rotation of 2M0036+18 then using 
R s ~ 6.3 x 10 9 cm, the inferred period indicates a surface equa- 
torial rotation velocity v w 37 km s -1 . This value is nearly 2.5 
times higher as that inferred from high-resolution optical spec- 
tra, vsin; w 15 ±5 km s" 1 (Schweitzer et al. 2001), suggesting 
an inclination angle, i w 24 ± 8°. Alternatively, if the period 
is related to the orbital motion of a companion which induces 
magnetic activity (see below) then the semi-major axis is about 
3. 1 x 10 10 cm or about five times the stellar radius; this is simi- 
lar to the case of the highly active RS CVn systems (e.g., Mutel 
& Lestrade 1985). Finally, the periodic rise in flux may be a 
series of weak flares, rather than variable persistent emission, 
in which case the flare generation process is periodic. 

7. THE ORIGIN OF THE MAGNETIC ACTIVITY IN 2M0036+18 

The results presented in the previous sections directly con- 
firm for the first time that the relative activity patterns in late-M 
and L dwarfs differ from those in early-type stars and early M 
dwarfs. It is therefore crucial to address the origin of the mag- 

19 We do not investigate time resolutions smaller than 2 min since the compul 



netic fields in the well-studied 2M0036+18 as an example of a 
larger trend that is emerging in sub-stellar objects. While some 
dynamo models for convective stars predict a reduced strength, 
physical scale, and lifetime of the magnetic field (e.g., Durney, 
De Young & Roxburgh 1993), the extent of this reduction is 
still not known. Clearly, the stability of the magnetic field of 
2M0036+18 over a period of three years indicates a long-lived 
process. Similarly, the physical scale of the field, R ~ R s , is sig- 
nificantly larger than the physical scale of the convection. Fi- 
nally, the inferred field strength and electron densities are sim- 
ilar to those inferred in early M dwarfs. It is therefore possible 
that the current turbulent dynamo theories are missing some key 
ingredients, or that the magnetic fields are generated and ampli- 
fied in another process, possibly by interaction with a close-in 
companion. 

7.1. The Influence of a Close-in Companion 

The enhancement of magnetic fields by a close companion 
may play a role in the RS CVn class of active short-period bi- 
naries, and a comparison is therefore particularly illustrative. 
These systems are known to produce radio gyrosynchrotron 
emission from magnetic fields of about 200 G, with a typical 
luminosity of L„. s w 2 x 10 16 erg s" 1 Hz" 1 , a factor of 10 3 
higher than that from 2M0036+18. Using the scaling for gy- 
rosynchrotron emission, L„ jS oc B~ 3 / 4 R 2 , and taking into ac- 
count the similar magnetic field strengths and covering frac- 
tions, as well as the ratio of radii (0.09 R vs. ~ 2 R Q ), we 
find L„,2M/£i/,RSCVn ~ 2 x 10~ 3 , in close agreement to the mea- 
sured ratio. In addition, the ratio L v ,r/L\, q \ for RS CVn sys- 
tems is known to be correlated with the orbital period of the 
system, such that tighter binaries are more radio active (Mu- 
tel & Lestrade 1985; Drake, Simon & Linsky 1989). Extrap- 
olating this relation to the value L v , R /L\, D i fa 6 x 10" 17 Hz" 1 
for 2M0036+18 we predict an orbital period of about 1 hour. 
This is not so dissimilar from the period of 3 hours inferred 
here. Thus, it is possible that 2M0036+18 represents a low- 
mass scaled version of the RS CVn binaries and that the same 
type of interaction that enhances the magnetic fields in RS CVn 
systems may take place here. 

Another illustrative example is the recent detection of en- 
hanced chromospheric Ca II H & K emission from the star HD 
179949 orbited by a 0.84 My planet with an orbital period of 
about 3.09 days, or a semi-major axis of 6.7 x 10 11 cm (Shkol- 
nik, Walker & Bohlender 2003). The emission appears to be 
periodic and in phase with the orbit of the planet. The conclu- 
sion drawn by Shkolnik, Walker & Bohlender (2003) is that the 
planet induces magnetic activity on the stellar surface which 
gives rise to chromospheric emission. 

These diverse examples illustrate that interaction with a com- 
panion may enhance activity. In this context both tidal and mag- 
netic effects can play a role (Cuntz, Saar & Musielak 2000). 
In the case of tidal interaction, if the orbital and rotation ve- 
locities are not the same, the generation of time-varying tidal 
bulges may give rise to increased turbulent motions or an en- 
hanced a-effect. The synchronization timescale (e.g., Goldre- 
ich & Soter 1966) for a companion orbiting 2M0036+18 is 
long, f sync « QMR 3 S /GM S )(M S /M p ) 2 (a/R s ) 6 w 1.6M s , y Gyr; a 
subscript p designates the primary, a subscript s designates the 
secondary, ui is the angular velocity, and a is the orbital separa- 
ting becomes prohibitively long. 
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tion. Thus, it is likely that the putative companion is not tidally 
locked and tidal interaction may be an effective mechanism. 

The strength of the tidal effect depends sensitively on 
both the binary separation and the mass ratio, Ag p /g p oc 
(M s /M p )[R p /(a-R p )] 3 . This suggests that the tidal effect may 
be particularly effective for late-M and L dwarfs where the pri- 
mary mass is low and the mass ratio may thus be of order 
unity. For example, for the HD 179949 system the distortion 
is Agp/gp « 1.3 x 10" 6 , while in the case of 2M0036+18 as- 
suming a Jupiter-mass companion and using a « 5R p we find 
Agp/gp « 3 x lO" 4 . 

Direct magnetic interaction is also possible, with the mag- 
netic energy of the system proportional to the product of the 
primary and secondary field strengths (Cuntz, Saar & Musielak 
2000). In this case the energy is E B oc B p B / 'd^pd^, where 
d mtS is the magnetospheric radius of the secondary and d m . p = 
a-d m . s is the distance from the primary to that radius. The en- 
ergy flux is related to Eb via F oc eEgVc, where e is an efficiency 
factor and v c is the combined velocity of the stellar magnetic 
motions (or turbulence) and the relative velocity of the orbit (at 
the stellar surface) compared to the rotation of the star. Thus, 
as in the case of tidal interaction, the possible larger mass ratio 
and small separation in the case of 2M0036+18 suggests that 
this mechanism may be more efficient than in extra-solar plan- 
ets orbiting solar- type stars. 

A final possibility is that a brown dwarf companion to 
2M0036+18 may undergo sustained Roche lobe overflow 
which may result in enhanced activity (Burgasser et al. 2000, 
2002). This scenario has been proposed for Ha active dwarfs, 
although no direct evidence was found to support this hy- 
pothesis. Here we may for the first time have the required 
evidence. For sustained overflow to occur, the logarithmic 
change in radius of the secondary, dlnRi/dlnM s < -1/3, re- 
quires a mass ratio q = M s /M p < 0.63. The period is given 
by P 2 = 4-TT 2 a 3 /GM p (l+q) while the radius of the secondary is 
equal to the Roche lobe radius, R s =R L = 0A9aq 2 / 3 /[0.6q 2/3 + 
ln^ + g 1 / 3 )]. The radius of the secondary is also given by R s = 
(4/3MqM p )- 1 / 3 /[l + 1.8(qM p T l/2 ] 4/3 (Zapolsky & Salpeter 
1969; Stevenson 1991), where the constant of 4/3 is appropri- 
ate for an age of 1 - 5 Gyr (Burgasser et al. 2002). Using 
P = 3 hr and M p = 0.076 M Q we find q « 0.14, consistent with 
sustained mass transfer, a secondary mass M s « 13 My, a sec- 
ondary radius R s « 1.1 Ry, and an orbital separation a w 5.2R p . 

Is there a way to directly test the hypothesis of a 
close-in companion? As in the case of planetary com- 
panions around solar-type stars, radial velocity, astromet- 
ric shift, and photometric techniques may provide an an- 
swer. The expected radial velocity amplitude (assuming 
zero eccentricity) in the case of 2M0036+18 is A w 2.3 x 
lO 3 (P/3hr)- 1 / 3 (M J sin/ o /M y )(M / ,/O.O76M ) 2 / 3 m s" 1 for a 
Jupiter-mass companion. Using a cross-correlation of the four 
strongest and narrowest absorption features of 2M0036+18 
(Figure 12) we obtain a limit of A < 4 km s" 1 from the pair 
of MIKE spectra separated by about half a period (Figure 13). 
This indicates that M s sini < 1 .7 My if the limit is on the max- 
imum amplitude. Naturally, if the two epochs were obtained 
when the putative companion was located ~ 90° away from the 
line of sight to 2M0036+18, the lack of a radial velocity signa- 
ture does not provide a useful limit. 

For an orbital inclination cos/„ < (R p + R s )/a, or i ;> 68°, 



we expect to observe an eclipse as the companion occults 
2M0036+18 (z' = 0° corresponds to a face-on orbit). For 
i a £ 73° more than half the surface of 2M0036+18 would be 
occulted resulting in large photometric variability. If all incli- 
nations are equally likely, the probability of eclipse is about 
24%. Gelino et al. (2002) find no significant variability for 
2M0036+18 in their I-band study of photometric variability in 
L dwarfs. Observations were obtained with a time separation 
as short as 1 - 2 hours, while the maximum time baseline in- 
vestigated was 53 days. The lack of variability indicates that 
i a < 68° //"2M0036+18 is in fact orbited by a close companion. 

Finally, a companion will induce an astrometric shift given 
by 9 « 0.lM s M 2 p /3 (P/3hr) 2 / 3 mas for a 2 : 1 mass ratio. This 
is just within the reach of current Very Long Baseline Interfer- 
ometry (VLBI) observations, but lower mass ratios will result 
in an undetectable shift. In principle, an astrometric shift down 
to a mass ratio of 50 : 1, or M s « 1 .7 My, may be detected with 
the 4 juas precision of the Space Interferometry Mission (SIM). 

7.2. Magnetic Structure and Rotation 

The 3 -hour period may alternatively be due to the rotation of 
2M0036+18. Using the inferred radius of 0.09 R , the equa- 
torial velocity is 37 km s" 1 , leading to an inclination angle of 
24 ± 8°, or nearly pole-on. The azimuthal orientation, <j>, of the 
axis of rotation relative to the line of sight is not known. How- 
ever, given the low inclination and the large covering fraction, 
the radio-emitting region has to be located at a low latitude for 
most values. This is because a location near the pole would 
make the bulk of the region visible throughout the rotation pe- 
riod, resulting in a nearly constant flux level. As can be seen 
from Figure 2 the flux drops to nearly zero between the peaks 
indicating the emission region is fully occulted for about half 
the rotation period. We note also that the low inclination may 
explain the constant sign of the circular polarization since for 
most values of <j> only one of the hemispheres is visible. 

If the 3-hour period is related to the rotation of 2M0036+18 
then we can estimate the Rossby number, Ro = P/t c , which 
is relevant for dynamo models. The convective turnover time 
for 2M0036+18 is r c = {MR 2 /L) l l 3 « 0.8 yr indicating that 
Ro i=a 4.4 x 10" 4 . At such low values early M dwarfs exhibit sat- 
urated X-ray emission, L x /L bo \ ~ 10" 3 (Pizzolato et al. 2003). 
Clearly, the conditions in L dwarfs are sufficiently different that 
a low Rossby number does not result in X-ray emission. We 
also note that the large value of t c may explain the stability of 
the radio emission (flux and circular polarization) over a period 
of about three years. 

7.3. Periodic Flaring 

Finally, it is possible that the variable persistent emission is in 
fact a series of periodic flares with a rate of occurrence of about 
0.33 hr" 1 . The frequency of strong flares, similar to the one 
detected by B01 with a flux of about 720 /zJy, is £ 0.04 hr" 1 . 
This pattern is similar to the case of the M9 dwarf LHS 2065 
for which strong Ha flares have an occurrence rate of <; 0.03 
hr" 1 while weak flares occur at a rate as high as about 0.5 
hr" 1 (Martin & Ardila 2001). If this is the case here, then the 
timescale to build up magnetic stresses is about 3 hours, while 
the lifetime of the electrons, corresponding to the width of the 
flares, is about 1 hour. 

In this scenario it is possible that the periodic weak flares 
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are inefficient at heating the corona and chromosphere, thus ex- 
plaining the lack of accompanying X-ray and Ha emission. The 
more rare strong flares, on the other hand, with an impulsive en- 
ergy release a few times larger may be accompanied by efficient 
heating. However, we note that the time-integrated energy re- 
lease in the putative flares observed here, E ~ vL v t ~4x 10 26 
erg, is not so different from the integrated energy release in- 
ferred for the stronger flare (B02). 

8. EMERGING ACTIVITY PATTERNS OF LATE-M AND L DWARFS 

The growing sample of late-M and L dwarfs observed in the 
radio and X-rays allows to investigate the patterns of activity 
as a function of relevant physical parameters. Ha observations 
of late-M and L dwarfs indicate a significant drop in the level 
of emission compared to the bolometric luminosity (Figure 14). 
Objects earlier than about M6 tend to have a ratio L Ha /L^ \ in 
the range of 10" 4 to 10~ 3 . However, by spectral type LO this 
ratio is typically lower than 10~ 5 . The few late-M and L dwarfs 
which exhibit flares (<~ 1%; Liebert et al. 2003) appear to have 
levels of emission similar to those of early M dwarfs, with only 
a single source exceeding that level to date. Thus, Ha emis- 
sion drops significantly with a decreased surface temperature 
(later spectral type), and flares exceeding the saturation value 
are extremely rare. 

Similarly, the ratio Lx/h, Q \ drops from a saturation value of 
about 10~ 2 5 in mid-M dwarfs to less than 10" 4 for the few late- 
M dwarfs observed to date (e.g., Fleming et al. 1993; Rut- 
ledge et al. 2000; Martin & Bouy 2002). With the exception 
of a recently-discovered M9 dwarf for which a ratio of about 
0.1 was measured during a flare (Hambaryan et al. 2004), 
the flares observed in other late-M dwarfs do not exceed the 
saturation level observed in the mid-M dwarfs (Stelzer 2004). 
Thus, in both Ha and X-rays the saturation effect appears to be 
real, with flares replacing persistent emission in late-M and L 
dwarfs. 2M0036+18 appears to follow the general trend with 
upper limits on the X-ray and Ha emission that are about three 
orders of magnitude below the respective saturation values. 

On the other hand, the radio activity in late-M and L dwarfs, 
and in particular 2M0036+18, exhibits a completely different 
trend. The observed ratio L R /L^\ is about 10~ 8 for persistent 
emission from early and mid-M dwarfs. However, the level of 
radio emission from late-M and L dwarfs exceeds this value 
by about an order of magnitude, while flares are brighter by at 
least two orders of magnitude. 2M0036+18 in particular has a 
ratio L R /L ho \ of about 10" 6 2 . Thus, while a significantly smaller 
number of objects have been observed in the radio compared to 
Ha, every single detection is brighter than the value observed 
in mid-M dwarfs. Moreover, the fraction of detected objects 
is much higher, <~ 40%. This suggests that a saturation effect 
does not exists in the radio band. Moreover, the radio luminos- 
ity does not seem to decrease with surface temperature. 

Another important trend, proposed by B02, is that rapid rota- 
tors exhibit stronger radio activity. This is shown in Figure 15. 
In Ha there is a clear rotation-activity relation in dwarfs earlier 
than M7, but this relation clearly breaks down in late-M and L 
dwarfs. A similar relation exists in X-rays, but it again breaks 
down in late-M and L dwarfs. In the radio band, on the other 
hand, the rotation-activity relation appears to extend to late-M 
and L dwarfs. The only objects for which the ratio of radio 
to bolometric luminosity is lower than the detected objects are 
those with vsin; <; 10 km s" 1 . 



Thus, it appears that the efficiency of magnetic field genera- 
tion and dissipation in fact does not decrease in late-M and L 
dwarfs, but the efficiency of chromospheric and coronal emis- 
sion does. Similarly, the rotation activity relation appears to still 
hold, but the activity is manifested in radio emission instead of 
Ha or X-rays. The detection of rare Ha and X-ray flares does 
suggest that chromospheric and coronal heating may still occur, 
but this may require large reconnection events, possibly accom- 
panied by large radio flares. 

9. CONCLUSIONS 

The simultaneous, multi-wavelength observations presented 
in this paper provide unparalleled insight into the magnetic field 
properties of an L dwarf. Most importantly, these observations 
directly confirm that radio activity is more prevalent in late-M 
and L dwarfs compared to Ha and X-ray activity. The detec- 
tion rate in the radio may be as high as 40%, while in Ha it is at 
most a few percent. Only a few objects later than M7 have been 
observed in the X-rays, but persistent emission is clearly rare. 
A detailed analysis of the radio emission along with the lack 
of detectable X-ray and Ha emission give rise to the following 
inferences: 

1. The radio emission from 2M0036+18, as compared to 
the bolometric luminosity, is the strongest observed 
in any dwarf star to date. Moreover, the emission 
from 2M0036+18 violates the radio/X-ray correlation 
by at least four orders of magnitude suggesting a dras- 
tic change in the coronal and chromospheric conditions 
compared to those observed in active stars down to spec- 
tral class M7. This is the first direct confirmation of this 
violation which has been observed in two additional late- 
M dwarfs (B01;B02). 

2. The radio emission is due to gyrosynchrotron radiation 
in a magnetic field of about 175 G. The size of the emis- 
sion region is about 10 10 cm (~ 1 AR S ), indicating a large 
covering fraction at the surface. These characteristics are 
similar to those inferred in early M dwarfs. 

3. The constant fraction and sense of the circular polariza- 
tion in observations conducted over a period of more 
than three years, implies that an ordered magnetic field 
is maintained stably on long timescales, despite the oc- 
casional production of strong flares (B02); the incidence 
rate for the latter is < 0.04 hr" 1 . 

4. The radio emission is periodic (^ 99.999% confidence 
level) with a period of 3 hours. If related to the rota- 
tion period this indicates a rotation velocity of about 37 
km s" 1 , while if related to an orbital period it indicates 
a semi-major axis of about five times the stellar radius. 
Alternatively, the period may be due to periodic release 
of the magnetic energy in the form of weak radio flares. 

Clearly, the radio emission from 2M0036+18 requires a 
large-scale, long-lived, and relatively strong magnetic field. 
The lack of accompanying Ha and X-ray emission suggests that 
the general decline in these activity indicators in late-M and L 
dwarfs is not a reliable tracer of decreased magnetic activity. In 
particular, if the magnetic field is essentially decoupled from 
gas in the cool atmospheres of late-M and L dwarfs so that 
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magnetic stresses are suppressed, this may quench the chro- 
mospheric and coronal emission. However, at larger distances 
the rate of collisions between charged and neutral particles may 
drop and the field could then couple effectively to the gas (Mo- 
hanty et al. 2002). This is the region where radio emission is 
expected to be generated most effectively, explaining the strong 
radio signal compared to the reduced Ha and X-ray emission. 
However, even this scenario cannot easily account for the severe 
violation of the radio/X-ray correlation, or for the detection of 
Ha and X-ray flares from some L and T dwarfs. 

While the origin of the observed periodicity remains un- 
known, all three scenarios provide interesting constraints for 
turbulent dynamo models. In the case of periodicity due to 
a close-in companion, it is possible that the magnetic field 
of 2M0036+18 is amplified by interaction with the compan- 
ion. This may alleviate the need to generate strong, long-lived 
fields from convection alone. If the incidence of activity in 
the radio band is truly 40%, this may also require an unusu- 
ally high fraction of companions, especially in comparison to 
current estimates of the binary fraction in late-M and early L 
dwarfs of ~ 15% (Close et al. 2003; Gizis et al. 2003). How- 
ever, this fraction is estimated primarily from direct detection 
with AO systems and the Hubble Space Telescope which are 
only sensitive to companions at separations of ;> 1 AU, sev- 
eral hundred times larger than the possible orbit in the case of 
2M0036+18. Direct evidence for close-in companions is re- 
quired before strong claims can be made, but radio active late- 
M and L dwarfs may provide a signpost for binary systems. 

On the other hand, in the case of periodicity due to rotation, 
it is clear that any dynamo model has to explain the presence 



of magnetic structures similar to those of M dwarfs without the 
benefit of a radiative-convective transition zone and without hy- 
drogen burning as a source of heat. Obviously, if the observed 
radio emission is in fact due to a series of episodic flares, this 
places stringent constraints on the timescale of magnetic field 
amplification and dissipation (~ 1 hour) over a large fraction of 
the stellar surface. 

The multi-wavelength approach taken in this paper is clearly 
the first step in addressing the details of magnetic field gener- 
ation at the bottom of the main sequence and beyond. Contin- 
ued studies, particularly in the radio and X-ray poorly-sampled 
range M9 to L5 (Figure 14), are essential for understanding in- 
dividual objects, as well as the effects of various physical prop- 
erties (e.g., binarity, rotation, temperature, age). These obser- 
vations will for the first time provide robust constraints for dy- 
namo models and for the structure of late-M and L dwarfs all 
the way from the interior to the corona. 
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APPENDIX 

X-RAY OBSERVATIONS OF THE L5 DWARF 2MASS 115074769-1627386 

In addition to the X-ray observations of 2M0036+18 outlined in this paper, we also observed the L5 dwarf 2MASS J15074769- 
1627386 (hereafter, 2M1507-16; Reid et al. 2000; Kirkpatrick et al. 2000) as part of our Chandra A03 program. This object is 
located at a distance of 7.4 pc (Reid et al. 2000, 2001; Dahn et al. 2002), has a surface temperature of about 1630 K, a bolometric 
luminosity of 10~ 427 L Q , an inferred radius of 0.09 Rq (Vrba et al. 2004), and a rotation velocity vsim w 27 km s" 1 (Bailer-Jones 
2004). The upper limit on Ha emission is L Ha /L\, \ < 10" 5 76 (Reid et al. 2000) and no Li is detected to a limit of 0.1A. Finally, 
B02 find an upper limit on the radio luminosity at 8.5 GHz of 3.8 x 10 12 erg s" 1 Hz" 1 from a 2.2-hr observation, corresponding to 
L R /L boX £ 10" 68 . 

The details of the X-ray observation and data analysis are given in §2.2. No X-ray source is detected at the position of 2M1507-16. 
The number of counts in the low and total energy bands were 2 and 4 respectively, while we expect 2 and 5 counts from background 
alone. We find no evidence for variability: In the total energy band, the shortest time-span which includes 3 counts is 2330 s, which 
we find from a Monte Carlo simulation, will occur in 22% of 27.7 ks light curves containing four counts. The 90% confidence limit 
on the time-averaged flux (<7 counts per 27.7 ks) is < 1 x 10~ 15 erg cm" 2 s" 1 . Taking the 3 counts within 2330 s, we place a 90% 
confidence upper limit on a 1-hr peak flaring flux of < 1 x 10" 14 erg cm" 2 s" 1 . At the distance of 2M1507-16, the upper limit on 
persistent emission is L x <; 6.6 x 10 24 erg s" 1 , or log(Lx/£*oi) <; -4.5. 

The upper limits on X-ray, radio, and Ha activity for 2M1507-16 are shown in Figures 14 and 15. Clearly this object follows 
the general trend of decreased Ha and X-ray activity in L dwarfs. The radio activity is at least a factor of four lower than that of 
2M0036+18, but the upper limit is several times brighter than other late-M and L dwarfs with a similar rotation velocity. 
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Fig. 1.— Radio light curves of 2MASS J00361617+1821 104 at 8.5 GHz (top) and 4.9 GHz (bottom). Both total intensity (black 
circles) and circular polarization (gray diamonds) are shown, with a 10-min time resolution. The negative sign of the circularly- 
polarized flux indicates left-handed polarization. The light curves for both frequencies and polarizations are variable, though with 
lower significance at 8.5 GHz (§6). 
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Fig. 2.— Fractional deviation relative to the mean of the circular polarization (bottom) and total intensity (middle) light curves of 
2MASS J00361617+1821104 at 4.9 GHz. Time resolutions of 10, 15, and 20 minutes are plotted. The top panel shows the light 
curves of two field sources with a 10-min time resolution. Clearly, the emission from 2M0036+18 is highly variable independent of 
time resolution. The light curve of 2M0036+18 also appears to be periodic with a period of about 3 hours. 
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Fig. 3.— Combined 4-hour optical spectrum of 2M0036+18. The location of Ha and prominent absorption features are marked. 
Clearly, no Ha emission is detected above the continuum level (inset). 
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Fig. 4.— Radio versus X-ray luminosity for stars exhibiting coronal activity. Data for late-M and L dwarfs are from B01 and B02, 
while other data are taken from Giidel (2002) and references therein. Data points for the Sun include impulsive and gradual flares, 
as well as microflares. The strong correlation between L« and Lx is evident and extends to spectral type M7 (see inset). Clearly, 
2MASS J00361617+1821 104 and the other late-M and L dwarfs detected in the radio to date violate the correlation by many orders 
of magnitude. 
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Fig. 5.— Auto-correlation function (ACF) as a function of temporal lag for the total intensity (black circles) and circular polarization 
(gray diamond) light curves of 2MASS J00361617+1821 104, as well as two field sources (gray squares). The time resolution is 10 
minutes (thick lines and field sources). A clear peak in the ACF of 2M0036+18 is seen at a lag of 3 hours. The flat ACF of the field 
sources is indicative of a non-periodic signal. The timescale of about 1 hour on which the signal from 2M0036+18 de-correlates (i.e., 
the first minimum in the ACF) indicates that the covering fraction of the emission region is ~ 10%. 
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Fig. 6. — Lomb-Scargle periodogram for the circularly polarized emission at 4.9 GHz with a time resolution of 8 min (black) and 
7 min (gray). The dashed lines indicate the 68.3%, 95.5%, 99.7%, 99.9%, and 99.99% significance levels based on Monte Carlo 
simulations of 50,000 randomized versions of the light curve (inset). For both time resolutions the most significant peak (99.9% 
and 99.99%) is at a frequency of 0.33 hr" 1 , or a period of 3.0 hours. We find the same result, though with different significance (see 
Figure 7), for both total intensity and circular polarization at a wide range of time resolutions. 
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Fig. 7.— Significance of the 3-hour period in the total intensity (black) and circular polarization (gray) 4.9 GHz light curves as a 
function of time resolution. The significance is determined using Monte Carlo simulated and randomized light curves as described 
in §6 and shown in Figure 6. 
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Fig. 8.— Follow-up observations of 2M0036+18 at 4.86 GHz taken on 2005, Jan. 10.99 UT for a total of 5 hours. The 3-hour period 
detected in the initial observations from Sep. 2002 is clearly seen. 
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Fig. 9.— Follow-up observations of 2M0036+18 at 4.86 GHz taken on 2005, Jan. 1 1 .99 UT for a total of 5 hours. The 3-hour period 
detected in the initial observations from Sep. 2002 is clearly seen. 
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Fig. 10.— Auto-correlation function (ACF) as a function of temporal lag for circular polarization in the follow-up observations on 
Jan. 10 (top) and Jan. 11 (bottom). The time resolution ranges from 2 to 15 min. A clear peak is seen at a lag of 3 hours, with a 
fall-off time of about 1 hour, confirming the periodicity. 
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Fig. 11.— Phased light curve (total intensity and circular polarization) for the observations from Jan. 10 and 1 1 UT folded with a 
period of 184 min. 
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Fig. 12. — Zoom-in on four spectral features in the high-resolution MIKE spectra of 2M0036+18, smoothed by two times the 
instrumental resolution. The second spectrum has been offset downward for clarity. No obvious shift in the lines is discernible. 
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Fig. 13. — Cross-correlation of the two MIKE spectra using the orders in which strong absorption lines were detected. The peak is at 
zero shift indicating no detectable radial velocity signature between the two epochs. The limiting factor in the accuracy of this result 
is the width of the absorption features. 
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Fig. 14. — Radio (top), X-ray (middle) and Ha (bottom) activity as a function of spectral type. Shown are objects with flares (solid 
squares) and persistent emission (detections: open squares; upper limits: triangles), as well as objects with persistent emission only 
(detections: circles; upper limits: triangles). The solid line in each panel marks the maximum level of activity for spectral types 
earlier than M7, and the dashed lines are an extrapolation to later spectral types. In both X-rays and Ha the level of persistent activity 
drops significantly beyond M7, while flares occasionally reach the same level of activity of early M dwarfs. In the radio band, on the 
other hand, several objects exhibit activity levels much stronger than that observed in early M dwarfs. This is particularly pronounced 
inthecaseof2M0036+18. 
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Fig. 15. — Radio (top), X-ray (middle) and Ha (bottom) activity as a function of rotation velocity. M dwarfs up to spectral type M8 
are shown in gray symbols, while dwarfs later than M8 are shown in black symbols. The early M dwarfs exhibit a clear rotation- 
activity relation in Ha and X-rays (gray dashed lines). This relation breaks down in the late-M and L dwarfs with many fast rotators 
exhibiting no Ha or X-ray emission. In the radio, on the other hand, faster rotation appears to be correlated with activity even for 
dwarfs later than M9. The large open circle designates the L5 dwarf 2M 1507-16 (see Appendix). 



